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Topological Majorana fermion (MF) quasiparticles have been recently shown to exist in semicon- 
ductor quantum wires with proximity induced superconductivity and a Zeeman field. Although the 
experimentally observed zero bias tunneling peak and a fractional ac-Josephson effect can be taken 
as necessary signatures of MPs, neither of them constitutes a sufficient "smoking gun" experiment. 
Since one pair of Majorana fermions share a single conventional fermionic degree of freedom, MPs 
are in a sense fractionalized excitations. Based on this fractionalization we propose a tunneling 
experiment that furnishes a nearly unique signature of end state MPs in semiconductor quantum 
wires. In particular, we show that a "teleportation"-like experiment is not enough to distinguish 
MPs from accidental conventional zero energy states, but our proposed tunneling experiment, in 
principle, can make this distinction. 
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Majorana fermions [I] are localized particle-like neu- 
tral zero energy states that occur at topological defects 
and boundaries in superconductors. Zero energy simply 
means that they occur at the chemical potential. A MF 
creation operator is a hermitian second quantized oper- 
ator = J which anti-commutes with other fermion 
operators. The hermiticity of MF operators implies that 
they can be construed as particles which are their own 
anti-particles [Tj-^. The key issues at this time in the 
condensed matter context are two fold, first, we must 
predict and characterize materials supporting Majorana 
fermions and second, we must detect them experimen- 
tally. In this paper we address the second issue of exper- 
imental detection by proposing a nearly sufhcent experi- 
mental signature for Majorana fermions. 

MFs have recently been proposed to exist in the topo- 
logically superconducting (TS) phase of a spin-orbit (SO) 
coupled semiconductor 2D thin film [U [5] or ID nanowire 
[5H7] with proximity induced s-wave superconductivity 
and Zeeman splitting from a sufficiently large magnetic 
field. In principle, the MFs in such systems may be 
detected either by measuring the zero-bias conductance 
peak (ZBCP) from tunneling electrons into the end MFs 
[3 [HI IS] or by detecting the predicted fractional ac 
Josephson effect [51 [71 [TUHT^ ■ The semiconductor Majo- 
rana wire structure, which will be the system of our focus, 
is of particular present interest since there is experimen- 
tal evidence for both the ZBCP [T3lfl6] and the fractional 
ac Josephson effect in the form of doubled Shapiro steps 

Despite their conceptual simplicity, neither the ZBCP 
nor the fractional ac-Josephson effect experiments con- 
stitute a sufficient proof of MFs at the ends of topologi- 
cal superconducting wires. A non-quantized (2e^//i) zero 
bias peak, such as that observed in the recent experi- 
ments |14H16| . can in principle arise even without end 
state MFs [TTHHI. Similarly, while the observation of 



doubled Shapiro steps are consistent with the TS state 
[2UH^ , a fractional ac-Josephson effect can exist even in 
Josephson junctions made of ordinary quasi-lD p-wave 
superconductors such as organic superconductors [llj or 
the non-topological phase of the semiconductor nanowire 
[24] . Given these caveats as well as the considerable com- 
plexity of existing experiments, there have been several 
alternative proposals to detect the presence of MFs [251 - 
[2S|. We emphasize that no single experiment can be 
expected to give a definite answer, so it is valuable to 
explore additional signatures that strengthen the case of 
MFs. Based on the inherent quantum non-locality of 
MFs [ini [23 HHHSI] : in this paper we propose an alter- 
native tunneling experiment on semiconductor Majorana 
wires that furnishes a nearly sufficient signature of end- 
state MFs. 

Non-locality arises in MFs because they differ from 
conventional complex (Dirac) fermions in that they have 
no occupation number associated with them. To define a 
quantum state of a system with MFs we must consider a 
pair of MFs. Take the pair of MFs -fa and 7^ at the ends a 
and 6 of a nanowire (NW) shown in Fig. [ll and construct 
zero-energy complex fermion operator (V — | [7a + i^b] 
associated with the pair of MFs 7a and 7b [lOj. The 
quantum state of the system is then determined by the 
eigenvalue of = d^d = 0, 1. Since is related to the 
MFs by 



we see that the state of the whole system is determined 
by non-local correlations between the fractionalized MFs 
7a and 7^. Our central concern is how to probe this 
non-locality to provide a robust and sufficient criterion 
of MFs. 

An immediate idea involves trying to inject an electron 
into 7a and retrieve it from 7f,. By connecting leads to 



the left and the right ends of the TS wire in Fig. [Tj one 
could imagine that an electron injected into the end a 
flips the occupation number Ud from = to = 1. 
The injected electron can then escape from the end b 
flipping the state back from Ud — \ to Ud = 0. Such a 
process where an electron can enter from one end a and 
exit at the other lead &, can be interpreted as a transfer of 
an electron, which we will refer to as Majorana-assisted 
electron tunneling. However, as has been discussed in 
previous works |291 132j , such a transfer occurs in a way 
so as to not violate locality and causality. 

The amplitude for the Majorana-assisted electron tun- 
neling [29l [30] can be written as 

S%'\T) = {9\la{r)i,m9), (2) 

where \g) is the initial state of the Majorana wire, r is 
the time-interval between the tunneling events and 7a, 6 
are the Majorana fermion operators at the left and right 
end of the wire. More generally the tunneling ampli- 
tude from the end a to the end 6, while the state of the 
nanowire system changes from g to g' , can be denoted by 

Sab'\'^)- When the state \g) is a ground state of the TS 
system with a deflnite fermion-parity Fp, the amplitude 
Slff^r) = —iFp is non-zero. In this case the amplitude 

S''^^\t) can be detected by a so-called coincidence mea- 
surement for the absence of an electron on the left lead 
and the presence of an electron on the right lead after 
the tunneling event [29] . 

However, a coincidence measurement does not directly 
imply a non-zero S^flf \T) in more general situations. 

The amplitude S^ff^ (r) in Eq. [ij reflects the amplitude 
for being able to transfer an electron from a to 6 while 
leaving the state \g) invariant. On the other hand, the 
measurement of the coincidence probability, Pc, does not 
keep track of the internal state of the system. For a 
general system (i.e. one that may be topological or non- 
topological) , Pc for an electron entering at a and exiting 
at b can be written more generally as 

Pc^J2 I(.92l7.7ll5i)l', (3) 

where 51, 32 £^re the internal states of the wire, which are 
not necessarily identical. Here for a non-topological sys- 
tem 7^ ^ are conventional fermionic operators. While TS 
systems with MFs have a non-degenerate ground state in 
a given fermion parity sector, more general systems with 
zero-energy Dirac end states may have multiple allowed 
values for 171 , 52 ■ Therefore, the coincidence probability 
Pc cannot be considered a unique signature for a topo- 
logical system. 

An important example of the inequivalence of Pc 
and S^^^^\t), to which we will return later, is a non- 
topological superconductor with a complex fermion zero 
mode at each end (such a system could arise because 
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FIG. 1. Schematic experiment to detect Majorana-assisted 
electron tunneling between the MFs a and b along the dashed 
lines. The electrons introduced in the system through the 
probe are transported via the Majorana-assisted tunneling 
process along the dashed line, which enclosed the flux "I>. The 
energy levels of these electrons, which oscillate as a function 
of flux are measured by the tunnel probe. By coupling 
the nanowires NW to the bulk superconductor SC through a 
Josephson junction J J one can introduce a gate-controllable 
charging energy proportional to the fermion-parity. 

we have an even number of weakly coupled Majorana 
fermions at each end). The quantum state is charac- 
terized by the occupancy na,ni, of the two conventional 
zero energy end modes. We can easily have Pc in 
this non-topological setup. Suppose the initial state is 
gi = {ua ~ = 0), then the sum for Pc in Eq. |3] would 
have a non-zero contribution from g2 = {ria — rih — 1). 
The tunneling of an electron from the lead into the zero- 
mode at a changes the occupation from ria = to = 1. 
On the other hand, the electron required to change the 
occupation of the state b from rih = to rih — 1 comes 
from breaking of a Cooper pair. The other electron from 
the broken Cooper pair is emitted into the lead in the 
vicinity near b. Note that the process conserves the num- 
ber of electrons within the system and cannot be elimi- 
nated even by the introduction of a finite charging energy 
|31j . Therefore in order to clearly distinguish this case 
from the process of Majorana-assisted electron tunneling 
(that also returns a non-zero Pc), we require S'ff\T) 
given m Eq. [2] itself to be non-zero. In other words, we 
require that the system return to the same state g after 
the tunneling process, so the same electron that enters 
at a leaves at b. In this paper we focus on how sjff"^ (t) 
can be directly measured to result in a unique signature 
of end state MFs [M] . 

As we now argue, the tunneling conductance experi- 
ment proposed in Fig. [l] only receives contributions from 
electron-tunneling processes where slffir) ^ and 

therefore can be used to measure 

^af\^)- The setup 
in Fig. [1] consists of an external quantum dot labeled 
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QD that is connected to the ends a and h through tunnel 
barriers. The gate potentials in the vicinity of QD are 
arranged so that a single confinement level at energy e 
is present near the fermi level /ip of the probe, which in 
turn is chosen to be close to the chemical potential of the 
electrons in the superconductor, /isc- The tunneling am- 
plitudes ta.b of the tunnel barriers connecting the state 
in QD to a, b respectively are chosen to be smaller than 
the energy separation between the energy levels of the 
QD, so that only a single QD level participates. This is 
possible even for substantially long quantum dots made 
of InAs because of the small effective masses of InAs. In 
the case where the nanowire system NW in the TS phase 
there are MF modes at the ends a and h. In this work 
we will restrict ourselves to the case where the system is 
well-described by the non-number conserving proximity- 
induced superconducting Hamiltonian, so that the effects 
of phase fluctuations, which might limit the velocity of 
the Majorana assisted tunneling process, are ignored. 

The physical picture is that the electron in the QD 
tunnels into the pair of MFs and effectively creates an 
interferometer such that the electron path encloses a flux 
$ through the loop. Thus our set-up effectively imple- 
ments a measurement of non-local correlation of MFs in 
nanowires, which is related to interferometric proposals 
to detect non-locality in MFs in topological insulators 
[5T] . However, the proposal in Fig. [I] is more suitable 
for implementation in semiconductor nanowire structures 
that have been realized in experiments. In the case that 
there is a finite tunneling amplitude to go from the left 
end a of the TS system to the right end, the energy level 
of the electron in the wire can be expected to develop 
an 2<I>o = hc/e periodic dependence on flux similar to 
the Aharonov-Bohm oscillations of the energy levels of an 
electron in a metallic ring. The probe, which is connected 
to the quantum dot QD by a weak tunneling matrix el- 
ement, is used to measure the shift in the quasiparticle 
energy level through a peak in the tunneling conductance 
between the probe and the superconducting lead SC . 

To demonstrate the connection between the electron- 
transfer amplitude across the ends a and b and the energy 
spectrum for electrons in Fig. [T] consider a Hamiltonian 
for the system labelled NW to be H^. To ensure that 
the analysis here applies to both the topological and non- 
topological cases for NW , we will label the low-energy 
manifold of states that are accessed by the experiment 
by \n). The TS system is coupled to the QD through a 
tunneling Hamiltonian written as 



V 



(tc.e^^%ct+/i.c). 



(4) 



where /q. are the fermionic states at the end of the TS sys- 
tem and is the creation operator of the active fermionic 
state in QD. The phases i^Jq are introduced by the vec- 
tor potential from the flux $. We will choose a gauge 
where ipy, = ^ and kpa = % = 27r<I>/2$o. In the absence 



of the tunneling perturbation V , the energy levels of the 
combined TS-QD system are written as Pn,no = Vn + noS 
where = 0, 1 is the occupation state of the fermionic 
state in the QD. Here i>n are the unperturbed energy lev- 
els of the TS and e is the energy of the fermionic level 
in QD. A subset Ug of the values taken by the TS state 
label n describes states of the TS system with energies 
near the ground state. 

The action of the fermion creation operators / takes 
the system from a ground state ng to a different state 
n. Assuming that the energy level e in the QD is not 
resonant with a fermionic level in the TS, the state n will 
have a different energy from the state Ug i.e. Pno,ng 
Pfto.ii- Integrating out these intermediate excited states 
then yields an effective Hamiltonian 



^no,n„,n' Prig .ng^rig ,n' 



E 



{no,ng\V\no,n){nQ,n'g\V\nQ,n) 

PfiQ.n 



(5) 



in the manifold of low-energy state ng,n'g with energy 
E 

~ Pna,n„- Note that for the perturbation theory to 
be valid it is necessary for the gate-tuned energy e to be 
tuned far off-resonance from quasiparticle excitation en- 
ergies in the TS wire by more than the tunneling matrix 
elements that depend on ta. 

The flux-dependence of the quasiparticle energy of QD 
measured in the ground state is written as 

Eg^{^) - E,p{0) = - i?„„=o,„,(<f) 

-i?„o = l,n.(0)+i?no=0,„,(0). (6) 

Using Eq.|4]the flux-dependent part of the energy for the 
uq = 1 QD state is determined from diagonalizing the 
Hamiltonian 



1.9."',) I I, |2o("s."s)/'_\ 



(7) 



where 5^'^'""'''' (t) = ('^g|//3(0)/i(r)|n^) is a generaliza- 
tion of the electron-tunneling amplitude defined in Eq. [2] 
to the case where the wire has multiple states. 

In the case where the NW system is topological with 
the operators fa,b being Majorana fermions, there is a 
unique ground state Ug in a given fermion parity sector 
FP. The Hamiltonian 5h in Eq. [7]is a scalar, so that the 
flux-dependent part of the energy can be written as 

5E{<^) ^ j dTe'^''-'^^[e'^'^s[%'>\T)+e-''^/^s'ff*{~T% 



(8) 

where (f) = 27r<I>/"I>o. Thus, the presence of a finite tun- 
neling amplitude S^^lf \T) defined in Eq.jijleads to a 2$o 
phase contribution to the energy 5E. 



4 




2*0 



FIG. 2. (a) Conductance as a function of probe voltage 
Vp and flxix $ in the presence of end MFs at a and b shows 
2$o = ^ periodic oscillations, (b) Conductance as a function 
of probe voltage Vp and flux $ in the absence of end MFs at 
a and b shows 4>o = periodic oscillations. 



In the case where the quantum state of the wire is as- 
sociated with the fermion parity Fp ^ and Fp^b of conven- 
tional zero-energy states at the left and right ends respec- 
tively, the quantum state is flipped (for example) from 
rig = {Fp^a = l,Fp^b - 1) to n'g = (Fp.b = ~l,Fp,a = 
— 1). As a result, only the off-diagonal in fermion parity 
matrix element terms in Eq. [7] are non-zero. This allows 
one to transform the (/)/2 phase-dependence to a </> phase- 
dependence by a gauge transformation. The resulting 
spectrum in the non-topological case is only $o-periodic. 
Thus the 2$o periodicity of the spectrum that may be 
observed for MFs is a signature of the fractionalization 
of the fermion parity operator Fp in a non-local way into 
fermionic operators defined at the ends [55] . 

The most direct way to measure the quasiparticle level 
E is by performing tunnel spectroscopy using a tunnel 
probe P (shown in Fig. [T]) as has been done in previous 
experiments. In the limit of weak coupling between the 
probe P the conductance measures the quasi-particle en- 
ergy level E and in principle the flux dependent part SE. 
This leads to a peak in the conductance at the quasipar- 
ticle energy that would be measured from the probe P 
and whose energy would depend on the applied flux with 
a periodicity of 2$o as seen in Fig. [2jb). In the absence 
of MFs, the superconductor can only transport electrons 



through Cooper pairs and the resulting level spectrum 
wiU be $o-periodic. 

However, the introduction of the tunnel probe P, intro- 
duces the possibility that the fermion-parity FP of the 
TS system can change by tunneling a fermion between 
the TS system and the probe. In fact, the fermion tunnel- 
ing process that would change the quasiparticle number 
is also the process that leads to a tunneling current that 
allows Eqp to be measured in this approach. By solving 
Eq. [?) one finds that changing FP in the topological su- 
perconducting case shifts the phase dependence 6E{^) 
by $0, so that allowing FP to flip randomly would re- 
turn the 2<I>o periodicity back to a $o periodicity and 
destroy the signature of the TS phase. The problem 
with FP flipping randomly can be fixed by introducing 
a charging energy Eq |31j . The charging energy makes 
the fermion-parity-changing process energetically unfa- 
vorable [3TJ[3S]. The charging energy term Ec = e^/2C, 
where C is the capacitance is introduced in the TS wire 
together with a Josephson coupling Ej ^ Eq to an infi- 
nite superconducting lead M, so that the superconduct- 
ing phase in the TS wire is essentially fixed apart from 
phase-slips [351 IMj • In this limit we find (see Supplemen- 
tary material for derivation) that the effective Hamilto- 
nian describing a pair of MFs is written as 



H, 



FpS + i{-fa{taCa + Kcl) + JbiUCb + tfccj,)), (9) 



where S 



Er 



(ic) 



3/4 



,-i^Ej/Ec and Fp 



(—1)^ = i^alb is the fermion parity operator. Here we 
are assuming that the wire is significantly longer than 
the decay length of the zero-energy end modes and at 
the same time short enough so that C is not infinite in 
which case 5 will vanish. 

In principle, it is possible to measure the 2$o period- 
icity of the quasiparticle energy indirectly without the 
probe P, which would change Fp. This can be done 
by measuring the periodicity of the total energy Etot in 
a fixed fermion parity sector Fp. While measuring the 
total energy Etot is difficult, it is possible to measure 
the current /($) = dq,Etot{^)- Such a measurement in 
Fig. [l] is identical to the fractional Josephson effect [TD] . 
As mentioned earlier, non-adiabaticity effects prevent the 
fractional Josephson effect from being an unambiguous 
signature for MFs |T . Fixing the fermion parity by using 
a finite charging energy allows us to use a DC measure- 
ment, in turn avoid the non-adiabaticity effects and thus 
to uniquely determine the topological character of the 
MFs. 

Summary and Conclusion: In this paper we have pro- 
posed a scheme for uniquely identifying the Majorana 
assisted non-local electron tunneling between two MFs 
at the ends of a wire in the TS phase. In principle, a 
non-local transfer of electrons in TS wires may be ob- 
servable by a coincidence measurement in the absence of 
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a charging energy of the TS system [551 130] ■ An experi- 
mental proposal for such a coincidence measurement has 
been given by using a joint measurement by two point 
contact detectors [37]. The non-local transfer of elec- 
tron in the presence of charging energy |31| , which fixes 
the fermion parity, can also be measured by a non-local 
conductance or transconductance between the leads , 
where a voltage applied to the left-lead a results in a 
current in the right lead b. Using results of Ref. for 
symmetric ta = t^, = t we find that 



dVa 



32K 



16r2 + (^2 _ Y2^2 + 8r2(52 + -1/2) ' 



(10) 



which clearly vanishes for (5 — > 0. Here F cx t2 is the lead- 
induced broadening of the MFs. However, as we have 
shown here that the Majorana assisted electron tunnel- 
ing process using either a coincidence detection [291 EZ] 
or by measuring the transconductance with a charging 
energy, while interesting, cannot be taken as a definitive 
signature of MF modes because even conventional near- 
zero energy states trapped near the spatially separated 
leads can also produce such non-local signature. 

In contrast to the possible measurements of non- 
locality discussed in the previous paragraph, our pro- 
posed non-local correlation experiment in terms of the 
flux dependence of the tunneling spectrum provides a di- 
rect verification of the non-locality of MFs in TS wires. 
We emphasize that the non-locality of the end state MFs 
arises from the non-locality of the fermion parity, which is 
unique to topological systems and cannot be emulated by 
conventional systems [33]. In our proposal fermion parity 
is fixed by a superconducting single-electron transistor 
configuration. The scheme proposed in this paper where 
the tunneling amplitude is measured as a flux dependence 
of an energy level should be suitable for testing for MFs 
in the semiconductor wire system-based experiments. 

Like the ZBCP and the fractional ac-Josephson signals, 
our interferometric signal will give a positive result if the 
system contains a single MF at each end. However, as we 
mentioned in the introduction, both the ZBCP and the 
ac-Josephson signal can be mimicked by non-topological 
physics. In our current proposal, considering the non- 
topological situation with two MFs at each end of a wire 
(as in Ref. the system can in principle be prepared 
in a non-local initial state so that our current proposal 
gives a positive result. However, as shown in the supple- 
mentary section, decoherence induced by coupling to the 
leads disfavors the non-local initial state and washes out 
any signature of non-locality in the present experiment. 
Hence, although a weakly coupled pair of MFs (or a con- 
ventional zero-energy state) at each end may still show a 
ZBCP at or near zero bias at finite temperature, our in- 
terference measurement will correctly demonstrate that 
the system is smoothly connected to a non-topological 
phase. 
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Derivation of the effective Majorana Hamiltonian 

To study the possibility of controlling the fermion par- 
ity variable n^, we consider a BCS Hamiltonian contain- 
ing a charging term on the TS wire that is characterized 
by a capacitance C. The TS wire is assumed to be cou- 
pled to an infinite superconducting lead M through a 
Josephson junction J. The resulting Hamiltonian for the 
system is written as 

H = J dxdx'ipl.{x)HQ{x(T, x'a')'i(ja-'{x') 

+ A[J dx^'^pl{x)^pl{x) - e-''^^{x)^^{x)] 

~ Ej cose + Ec{n + nF/2-ngf, (11) 

where Hq is the single-particle part of the Hamiltonian 
representing the semiconducting wire, A is the proximity- 
induced pairing potential, Ej is the Josephson cou- 
pling between the superconductor and the wire, Ec{n + 
np /'i — TigY is the charging energy resulting from a global 
Coulomb interaction within the superconducting-island- 
nanowire system. Here n is the number of Cooper pair 
on the island, which satisfies the canonical commutation 
relation [n, 9] = i and ng is the gate charge of the TS 
system that is controlled by the total gate voltage on 



the system, np = J dxip^ {x)'4>{x) is the total number of 
fermions on the wire. 

To solve the low-energy eigenstates of this problem, 
we first extend the Hilbert space for the variable 9 to 
allow the wave-function to be 47r periodic. Later we will 
reject all solutions in the extended Hilbert space that 
do not satisfy the 27r periodic boundary conditions. In 
the extended Hilbert space, it is convenient to apply a 
unitary transformation U{9) = e*"^^/^, which transforms 
^ U{9yip'^U{9) = V^e**/2_ Tj^g transformation U{9) 
also transforms the term n = —ide — > n — np/2. The 
transformed Hamiltonian is now written as 

H = j dxdx''ipl^{x)Ho{xa, x'a')ipa'ix') 
+ A[J dx^{x)ijl{x) - 

- Ejcos9 + Ec{n-ngf. (12) 
The Hamiltonian H separates into an Hbcs written as 

Hbcs = J dxdx'^l{x)Ho{xa, x'a')i)„'{x') 

+ A[J dx^plix)^lix) - ^p^ix)^^ix)], (13) 

which commutes with n, 9, together with a Hjj written 
as 

Hj = -Ejcos9 + Ec{n-ngf, (14) 

where [Hj, Hbcs] = 0. 

In addition to Hj and Hbcs it is necessary to consider 
the fermion tunneling terms from the leads which in the 
transformed basis is written as 

i,../..W.K"VWc(x, + M, (15) 

where c'{x) is a fermion creation operator in the lead. 
Specifically, we will be interested in the case where the 
BCS Hamiltonian has low-lying (potentially Majorana) 
states fl fj at the ends a and b. Eliminating the low- 
energy gapped modes one can write 

Ht^ J2 [e"'/'(ijv,a/I + tAMca. + h.c]. (16) 

a—a,b 

To return to the original Hilbert space where 9 is 27r peri- 
odic, it is necessary to apply the unitary transformation 
Ui{9) = e-"''^/^, where = Ea flU In the case 
where fa are Majorana fermions, one must define a con- 
ventional fermion d^ = fa + ifb and define n^; — d^d. In 
either case, the applied unitary transformation replaces 
the e'^/^ dependence of Ht by an e'^ dependence. The 
Hamiltonian for the system now becomes 

H -Ejcos9 + Ec{n + nd/2-ng)^. (17) 
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The spectrum of the above Hamiltonian can be gauge 
transformed so as to shift n by an arbitrary integer. 
Therefore, the spectrum of H is integer periodic. The 
ground state shift induced by is given by 5 w 

-^Ec (j^) g-5-\/-Ej/Sc^ which is derived by solving 
the Mathieu equation. 

Assuming that to,& ^ \/E!jEc, which is the excitation 
gap of the phase variable or equivalently the Josephson 
plasma frequency, one can project into the ground state 
of the phase variable. Since, Ej/Eq 1, the phase 
variable 9 is localized near ~ and under the projection 
to the ground state the Franck-Condon factor e'* can be 
approximated to be unity. 



Dephasing of non-local correlations/entanglement 
for conventional fermionic zero-modes 



The non-local signatures of the MFs discussed previ- 
ously were consequences of robust non-local correlations 
that were unique to MFs. In the case of conventional 
low-energy (even zero-energy) modes such correlations 
cannot be present in a robust way. In particular, even 
if the system accidentally reaches a non-local correlated 
state, coupling to either Bosonic or fermionic dissipative 
bath will destroy such correlations. In the present section 
of the appendix, we show that coupling to a Fermionic 
bath leads to such desctruction of non-local correlations. 
We accomplish this in two sub-sections. In the first sub- 
section, we show that when placed in contact with a 
fermion bath in thermal equilibrium, the system reaches 
a unique thermal equilibrium state. In the second sub- 
section, we show that the system in thermal equilibrium 
cannot have any non-local correlations of the type that 
might lead to the measured interference or non-local tun- 
neling amplitude. 



Equilibration in the non-interacting case 

The non-interacting fermion problem can be solved in 
terms of 



Pab{t) = i{-1a{t)lb{t)). 



(18) 



The Hamiltonian in terms of Majorana operators is writ- 
ten as 



a,b 

The equation of motion is written as 

daa = ha,blb, 



(19) 



(20) 



with ha^b being anti-symmetric. The time-evolution of 
7a (t) is written as 



7(i) = G(«'(i)7(0), 



where G'(^)(i) = e^* for t > 0. Expanding G^^\t) in 
eigenstates 

G(^)(t)=^VnV'ie-"* 



-<■ 



iS 



— ih — iS 



(22) 



The time evolution of the density matrix is written as 

p(t) = G(^)(t)p(0)G(-^)(t). (23) 

Assuming that we start with a density matrix that is 
diagonal between system (s) and bath (b) (i.e. lead), the 
system density matrix 

pM = Gi''Ht)pMGi^Ht) + Gi"b\t)PbiO)Gif\t). (24) 
Assuming the system-bath structure for the Hamilto- 



nian and defining 



iR,A) 



(uj^iS — Hb) ^ , one can write 



the final system Green function as 

G(^-^)(c.) = {lot^S- H,,gi'''^\Lo)H,,)-\ (25) 

The system bath interaction Green function can be writ- 
ten as 

Gif^)(^)=5f^'(a.)i/,.G(^^^^)(..). (26) 

Provided the anti-hermitean part of the fermion-self- 
energy 



(27) 



does not have a null-space (i.e. a space of zero- 
eigenvalues) the Green function G')i^''^\t) can be as- 
sumed to decay exponentially in time. Therefore it fol- 
lows from Eq. [24] that the first term must vanish. To 
evaluate the second term in Eq. [241 we note that from 



(28) 



Eq. 26 it is clear that 

Gif^^(i) = J dTGs{t-T)H,b9b{r). 

Since Gs(i — r) is exponentially decaying, the 0{t) com- 
ponent of gii is not relevant and one can approximate 



G 



-hbT 



/•oo 

/ dTGs{T)Hsbe~ 
Jo 



■hbit- 



(29) 



up to exponential factors. Applying this identity to 
Eq. [24|we obtain 



Ps{t) 



(21) 



dTirfT2Gf)(ri)iJ,fce-''''(*-^iVfc(0)e^'(*-"^)i?„sGi^)(T2) 

dTidr2G(«'(Ti)ff.6P6(0)e''^(^^-^=)i?b,G(^)(r2), 

(30) 
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which is a t-independent asymptotic value. Transforming 
to frequency space, 

/oo 
dujfiuj)Gi''\co)H,,A,{ij)Ht,,Gi^\uj). 
-OO 

(31) 

Considering the integrand at ui away from a pole of gs{ui), 

= zGi«)(c.)i7,4.gf '(a.) - gi^'^u^nm^Gi^^u;) 

= *G(«)(c.){E(^)(^)-S(^)(c.)}G(^)(c.) 

= *{Gf) (a.) -Gi^) (..)}, (32) 

which can be checked using Dyson's equation. Therefore, 
in the absence of a protected null-space in the dissipative 
part of the self-energy "Eg^uj), the density matrix equili- 
brates to the grand-canonical thermal equilibrium value 

/oo 
du;/(a;){G(^)(u;)-G(^)(u;)}. (33) 
-oo 

Correlations in the grand-canonical thermal state 

First we prove a general lemma: Consider the grand- 
canonical thermal partition function of a system of Ma- 
jorana fermions (or fcrmions) which can be partitioned 
into two parts L and R, so that the Hamiltonian has a 
symmetry ^L,a — > —jL,a for all MFs ■^L,a on the left 
half of the system L. Then all correlation functions 



{lL,alR,b) = 0. To prove this simply apply the trans- 
formation 7L^a — >• —UjL,aU^ tO 

{lL,alR,b) OC Tr[7L,a7fl,6e-^^] = Tr[UjL,aUhR,bUe-^''U^ 
= -Tr[jL,alR,be-^"] = 0. (34) 

Therefore all such non-local fermionic correlations must 
vanish in the thermal state. 

As discussed before, in the topological state, the charg- 
ing energy gives rise to a fermion parity dependent term 
in the Hamiltonian 

Hni oc i^alb, (35) 

which violates the conditions for the above lemma and 

leaves a non-zero value for this non-local correlator. 
Therefore the presence of the Majorana fermion term 
is crucial for the appearance of a non-local correlation 
function that can contribute to the interference. Such an 
interference will not happen even if the conventional end 
zero mode is made of a pair of Majorana fermions. This 
is because, as discussed in the previous sub-section, cou- 
pling Majorana fcrmions to leads causes them to thermal- 
ize into the grand-canonical ensemble. Once they have 
thermalized in this ensemble, Hni violates the conditions 
of the lemma and generates a non-local fermion corre- 
lation. The long-range Coulomb interaction for conven- 
tional fermionic modes docs not lead to the violation of 
the conditions of the lemma and does not lead to any 
long range correlation. 



